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The covariant technique for calculating the heat kernel asymptotic expansion for an 
elliptic differential second order operator is generalized to manifolds with boundary The 
first boundary coefficients of the asymptotic expansion which are proportional to t 1 ! 2 and 
t 3 / 2 are calculated. Our results coincide with completely independent results of previous 
authors. 



1. Introduction 



One of the most fruitful approaches in quantum field theory, especially in gauge the- 
ories and quantum gravity, is the approach of the effective action. It is applicable for 
topologically nontrivial manifolds (noncomplete manifolds with boundary etc.) as well. 
The effective action is calculable only within the limits of some perturbation theory. In 
general case of arbitrary background it is impossible to get an exact answer even in one- 
loop approximation. That is why one should use various approximate schemes for the 
investigation of various aspects of the effective action. One of the most important features 
of such calculational schemes should be manifest covariance of the calculations at each 
order. 

The general framework of covariant methods for calculating the effective action is the 
heat kernel method [1]. Heat kernel plays very important role in quantum field theory. It 
has been successfully applied to the analysis of the structure of ultraviolet divergences and 
anomalies and renormalization [2] as well as to calculation of the finite part of the effective 
action that is expressible finitely in terms of the coefficients of the heat kernel asymptotic 
expansion. 

Various methods for calculating heat kernel asymptotic expansion were developed [3- 
5]. However, most of them come to nothing more than manifolds without boundaries and 
can not be applied immediately to boundary problems. 

Recently some papers dealing with boundaries appeared. These are, mainly, indirect 
methods that use the factorization properties of heat kernel and its behaviour under con- 
formal transformations [6-8] or an expansion of the boundary in the neighbourhood of the 
tangent plane in some point [9,10]. On the other hand the general methods used in math- 
ematical literature [11] are applicable in any case but the general covariance is lost and 
the method becomes not effective at higher orders. Nevertheless, the most complete sum- 
mary of formulae concerning boundary contribution in heat kernel asymptotic expansion 
is presented in [6]. 

In this paper we propose a new algorithm for calculating boundary contributions 
in heat kernel asymptotic expansion. We show how the ordinary technique [1,4] can be 
applied for manifolds with boundary. A very close approach was developed in [12]. 



Let M be a <i-dimensional compact riemannian manifold with smooth boundary dM, 
and F be elliptic second order differential operator on M of the form 



where V is a connection on smooth vector bundle V over M and Q is an endomorphism 
of this bundle. 

The heat kernel is defined by the equation 



2. The general framework of calculations 



F = - □ +Q 



(1) 




(2) 
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with initial condition 

U(0\x,x') = g- 1/2 (x)5(x,x') (3) 

In the case of manifolds with boundary one should also impose suitable boundary 
conditions: 



BU(t\x,x') 

Dirichlet 



= (4) 

xedM 



5 = 1 (5) 

or Neumann ones 

B = (n^ + S) (6) 

where n /i (x) is the unit inward pointing normal to the boundary and S is an endomorphism 
of the bundle V. 

In the case of compact manifolds without boundary conditions of periodicity are used 
instead of the boundary conditions (4) . One can show that for any t > the heat kernel is 
a smooth analytic function on the manifold M, that behaves like a distribution (5-function) 
near the boundary. 

In this paper we shall investigate the heat kernel, mainly, when the points x and x' 
are close together, as it is the limit x — > x' which is of interest in quantum field theory for 
calculating the effective action and vacuum expectation values of local observables. 

Having in mind the calculation of heat kernel asymptotic expansion at t — > let us 
take advantage of the quasi-classical approximation. Let us represent the solution of the 
equation (2) in the form 

U(t\x, x') = (Ant)- d / 2 A 1 / 2 (x, x') exp (- ^V^ ) to(t\x, x') (7) 
where a(x, x') is a symmetric biscalar satisfying the equation 

^=2^=2^'^' . ( 8 ) 

and another biscalar A(x,x') is defined by 

A(x, x') = g- 1/2 (x) det (-V„ V u a(x, x')) g- 1/2 (x') (9) 
and satisfies, as it is easy to show from (8), the equation 

A -i/2 DA i/2 = I( rf _ D(7 ) (10) 

Making use of (2) and (7)-(10) we obtain a transfer equation for the function Q(t\x, x') 

(J^ + ^D + A- 1 / 2 FA 1 / 2 ^n(t\x,x')=0 , (11) 



D = a^V 



The initial and boundary conditions are to be determined from (3) and (5)-(6). But first 
let us discuss the general structure of this construction. As is generally known [1], the 
equation (8) determines the geodetic interval defined as one half the square of the length 
of the geodesic connecting the points x and x' . However in general case of topologically 
nontrivial manifolds there are more than one geodesic between points x and x', i.e. the 
equation (8) has more than one solution. Therefore the quasi-classical ansatz, in general 
case, should have a form of a sum of analogous contributions from all geodesies. There 
is always one leading solution a(x, x') which is determined by the shortest geodesic. It is 
marked out by the fact that it goes to zero 



when the points x and x' approach each other. By the square brackets we denote here and 
below the coincidence limits of two-point quantities when the points x and x' tend to each 
other along the shortest geodesic. 

If one fixes a sufficiently small region including the points x and x' (when they are close 
enough to each other) then only this single solution is left. Therefore it is in some sense 
local and does not depend on the global structure of the manifold. Multiple geodesies are 
closely associated with two reasons reflecting the essentially global (topological) aspects 
of the manifold. First, there could be manifolds with closed geodesies. In this case in 
addition to the shortest geodesic there are always geodesies that emanate from point x' , 
pass through the whole manifold one or several times and return to the point x. Second, 
geodesies could be reflected from boundaries of the manifold one or more times. 

In general case there could be infinitely large number of additional geodesies. They can 
be ordered according to the value of the geodetic interval. It is obvious that the more times 
the geodesic is reflected from boundaries or passes through the whole manifold the larger 
the geodetic interval is. And according to (7) the value of the geodetic interval immediately 
determines the weight of the contribution of each geodesic in quasi-classical approximation. 
Let us mention that in heat kernel asymptotic expansion contribute only those geodesies 
for which the geodetic interval could vanish. (Note that in the euclidean signature used in 
this paper all the geodetic intervals are non-negative according to definition.) There is only 
one such a geodesic among all additional ones, namely, the geodesic with one reflection 
from the boundary. It is evident that it has the minimal geodetic interval in comparison 
with all others. Moreover, if one fixes sufficiently narrow strip of the manifold near the 
boundary including points x and x' then the geodesic with one reflection will be the only 
additional geodesic. In that sense the corresponding solution to the equation (8) is also 
local, i.e. it reflects the local properties of the boundary and does not depend on the global 
structure of the manifold. The corresponding geodetic interval <f>(x, x') in the coincidence 
limit x = x' equals doubled square of the normal distance to the boundary and, therefore, 
vanishes on the boundary 




x—>x' 



(12) 




= 



(13) 



dM 
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Therefore the corresponding solution contributes in the heat kernel asymptotic expansion 
but only on the boundary. All other solutions to the equation (8) are essentially global 
and corresponding geodetic intervals are strictly positive and do not vanish anywhere. So 
for the analysis of the heat kernel asymptotic expansion it is sufficient to restrict oneself 
to the local term and the term with one reflection 



U(t) = (4,rt)-*/a (exp (~) A 1/2 0(t) + exp (-|) * 



(14) 



where is the corresponding transfer function. (Here it is convenient not to single out 
the pre-exponential factor, i.e the Van Vleck-Morette determinant.) The functions and 
\1/ satisfy equations analogous to (8)-(ll) 



1 



1 



L^> = 



(15) 



(16) 



Since the second term in (14) does not contribute in the limit t — > outside of the boundary 
and remembering that [A] = 1 we get from (3) the initial condition for O. 



[n(o\x,x')] = i 

The boundary conditions (4) take now the form: 
Dirichlet 



xedM 



Neumann 



(n^ + S) (a^q + *) - L a , n (A 1 /^ - * 



xedM 



(17) 



(18) 



(19) 



xedM 



where a n = n^V^cr. 

For calculating the function O it is sufficient to take advantage of the standard 
Schwinger-De Witt expansion 



m = E 



-t) k 



k=0 



(20) 



where a,k(x,x') are the so called Hadamard - Minakshisundaram - De Witt - Seeley coeffi- 
cients (HMDS) which are determined completely independently of the boundary conditions 
by recursion relations 

(21) 



(l + \d) a k = A- 1 /2 FA V2 afe _ 1 



5 



The initial conditions for these recursion relations are obtained from (17) and have the 
form 

[oq] = 1 (22) 

To solve these relations we have elaborated in the papers [4] a special technique. It allows 
to calculate an arbitrary HMDS-coefficient in a sufficiently effective way. It is shown there 
that the formal solution of the recursion relations (21) 



( 1 +l D ^j 1 A-^FA 1 /* + _1_£^ A- 1 /2 F A 1 /2...(i + jD )-i A -i/2 FA i/2 

(23) 

takes the practical meaning in the form of covariant Taylor series 



a * = V E { -^~ ^ ---^ [ V (M! • • • V ^)«*] ( 24 ) 

n>0 

where V(x, x') is the parallel displacement operator along the shortest geodesic from point 
x' to the point x. In [4] the coefficients of that series are calculated and the results of 
calculations up to are listed. That is why we will concentrate our attention here on the 
calculation of the function 

Let us choose for convenience of further calculations a special coordinate system in 
the neighbourhood of the boundary x^ = (r, l ), where r is the length of the geodesic arc 
normal to the boundary in the point 9, the equation of the boundary having the form 
r = 0, and 9 % are the normal riemannian coordinates on the boundary. 

We assume all quantities to be analytic in the coincidence limit on the boundary 
6 = 9'. Hence we will expand all quantities in covariant Taylor series in the neighbourhood 
of the boundary. 

The further strategy is rather simple. One should introduce a small parameter reflect- 
ing the fact that parameter t is small and the points x and x' are close to each other and 
to the boundary 

t 1 / 2 ~ r ~ r' ~ (9 - 9') ~ e 

and construct a corresponding perturbation theory in this parameter. We expand the 
transfer operator L (16) in a formal series in the small parameter e 

L (e 2 t\er, er' , e(9 - 9'), 9') = -^L_ 2 + ^L_ x + L + ■ ■ ■ (25) 
and seek for a solution to this equation of the form 

* (e 2 t\er, er',e(9 - 9'), 9') = (t\r, r', (9 - 9'), 9') (26) 

n>0 

From the transfer equation we get easy recursion differential relations 

L_ 2 ^o = (27a) 
L_ 2 ^! = -L_!^ (276) 

L_ 2 #2 = - L *o ( 27c ) 



n 
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etc.. 

One shows easily that the coefficients \l/ n satisfy the scaling properties 

* n (e 2 t\er, er', e(9 - 9'), 9') = e n V n (t\r, r', (9 - 9'), 9') (28) 

Analogous equations take place also for the function O. Whereas these relations are suf- 
ficient for the calculation of all fi n , for single-valued calculation of \l/ n it is necessary to 
use additionally the boundary conditions (18) (or (19)). 

The main difference between them is that fi n are analytic in all variables while \E f n 
are analytic in (9 — 9') but are complicated functions of the variables t, r and r' (one can 
show that they are analytic in the variables R = (r + r') and u = (r — r') / (r + r'), so that 
the point r = r' = is singular). 

Making use of the scaling property (28) one can get from (26) an important represen- 
tation for #(£) 



*(t\r,r>,(e-e>),e') = J2t n/2 *n U 

n>0 ^ 



r r> (9-9') / 

TTTF^r' > (29) 



Using (11) and (14) we write down for the trace of the heat kernel 
Tr[fU(t)}= f dxg 1 / 2 tr[fU(t)] = (30) 

J M 

= (4nt)- d / 2 tr I [ dxg^fm} + [ d9^/ 2 Trfrexp ^ [/*(*)] 

Um JdM JO \ 1 J 7 ' 

where f(x) is an arbitrary smooth function on the manifold and 7 = det7i 7 -,7y is the 
induced metric on the boundary. One should mention that additional terms due to the 
presence of the boundary contribute to the asymptotic expansion of the trace of the heat 
kernel at t — > only when integrating over a sufficiently narrow neighbourhood of the 
boundary. Therefore, the upper limit of the integration over the distance to the boundary 
r in (30) is taken to be 00 (up to exponentially small terms). 

At last, using (20) and (29) we obtain the asymptotic expansion of the trace of the 
heat kernel both in volume part and in surface one) 

/ dxg 1 / 2 J2 L n-[fa k ]+ d»7 1/2 £f^c^(/)hl) 
fc=o k - JdM k=o 2 J 

where 

/•oo r(n-m) 

c±±i(f)= £ / rf£exp(-£ 2 )£ n ^-n(£) J2 —,r v ( 32a ) 



b k (i) = t- k/2 [** 

/( fc )=n^---n^V Ml ---V Mfe /| r=0 



(326) 

r=£,Vt 



7 



d_ 

9m~ I 



,1/2 



7 



1/2 



(33) 



r=0 



Thus for the calculation of the boundary contributions in heat kernel asymptotic expansion 
one suffices to compute the coincidence limits [\P n ]> put t = 1, integrate with the weight 
exp(— r 2 )r m and combine them with the quantities gu (33) calculated before. 



3. Explicit expressions 

Here we list the results of calculations of the lower order coefficients of the heat kernel 
asymptotic expansion omitting most of the cumbersome computations. The most complete 
list of volume contributions is presented in [4] . The simplest first of them have the form 



[a ] = 1 



f 



[oi] = Q--R 

O 



(34) 



H = (q-I r ) ~ I n Q ~ ^ R ^ Rflv + ^ R ^R^ aP + ^ □ r + 

where TZ^ = [V M , V n ]. 

Using the set of normal coordinates introduced above it is not difficult to calculate 
the coefficients gj- (33) 



90 = 1 

91 = -K 

g 2 = -K ij K" +K 2 -R° m 



(35) 



where 



K = ~f ij K, 



3 2 dr 3 



r=0 



lij — 9ij 

r=0 

RO = n ^ n -R 



r=0 



For simplicity we list below the solution only for Dirichlet boundary conditions. From 
equation (27) and boundary condition (18) we have 



#n = -1 



(36) 
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The calculation of next orders is considerably more difficult though it offers no particular 
problems. The result has the form 



[* 2 ]=J(Q- 1 -R)- 1 -(l+ r 4)K 



6 / 3 V t 



+ h {^) K2+h i^) K " K "} (37) 

where all tensor quantities are calculated on the boundary and R is the scalar curvature 
of the boundary 



h( z ) = / dxexp(— x — 2zx) = ex.p(z )Erfc(z) 
Jo 

fi(z) = ± + ^(2+±z 2 -z(z 2 + 6)h(z 



(38) 



Using these expressions and the integrals 



Jo 



r ( n + 2 \ 

2\tn U (t\ V 2 J 



/ ^ * ^ vs/ 2(n+l) 
we get finally several first boundary coefficients in asymptotic expansion (31) 

<*/,</) = f { {-\& - \< n + |a' 2 - ^ + Q) f 



(39) 



These results coincide with ones obtained by completely independent methods in [6] that 
confirms that the approach developed in this paper is correct. 

Obtained results may be used when investigating Green function and the energy- 
momentum tensor of quantum fields near the boundary. 



3. Conclusion 
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In this paper a new method for calculating the heat kernel asymptotic expansion (or 
one-loop effective action in any model of field theory) in the case of arbitrary background 
fields and manifolds of nontrivial topology (with boundary) is proposed. This method 
is a generalization of our previously elaborated covariant technique for calculation of the 
effective action and is based on the summation in the quasi-classical approximation the 
contributions of all geodesies connecting any two points of the manifold (including the 
geodesies reflected from the boundary). It is established that taking into account one 
additional geodesic with single reflection correctly reproduces all boundary contributions 
in the heat kernel asymptotic expansion. Proposed approach allows not only to calculate 
the asymptotic expansion of the trace of the heat kernel but also to analyze the local 
structure of the heat kernel near the boundary. Furhtermore we are going to calculate by 
means of it also the next terms of the asymptotic expansion (31) ~ t 1 and ~ t 5 / 2 , i.e. the 
coefficients C2 and C5/2 (32). 
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